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ABSTRACT 

We test the hypothesis that the polarization vectors of flat-spectrum radio sources (FSRS) in 
the JVAS/CLASS 8.4-GHz surveys are randomly oriented on the sky. The sample with robust 
polarization measurements is made of 4155 objects and redshift information is known for 
1531 of them. We performed two statistical analyses: one in two dimensions and the other 
in three dimensions when distance is available. We find significant large-scale alignments 
of polarization vectors for samples containing only quasars (QSO) among the varieties of 
FSRS’s. While these correlations prove difficult to explain either by a physical effect or by 
biases in the dataset, the fact that the QSO’s which have significantly aligned polarization 
vectors are found in regions of the sky where optical polarization alignments were previously 
found is striking. 

Key words: polarization - galaxies:active - quasars: general - radio continuum: general - 
large-scale structure of Universe 


1 INTRODUCTION 


Hutsemekers et al. j 199811200 iL l2005h have reported intriguing 
alignments of polarization vectors of quasars at optical wave¬ 
lengths. The presence of these alignments may possibly be in ten¬ 
sion with the well-accepted concordance model of cosmology as 
they involve correlations between sources separated by comoving 
distances at the Gpc-scale. 

Paying p articu lar attention to instrumental biases, 
I Jackson et al] l l2007h compiled the JVAS/CLASS 8.4-GHz 
sample of flat-spectrum radio sources (FSRS) with polarization 
position angle (PA) measurements. As they have shown, rotation 
measures of the polarization vectors induced by Faraday rotation 
at 8.4 GHz are too small to destroy information about the intrinsic 
PA’s. Therefore, any observed correlation of PA’s among sources 
can be thought to be intrinsic to the sources themselves. Extracting 
from this sample 4290 FSRS’s with polarized flux higher than 1 
mjv. ljoshi et alJ ( l2007h searched for systematic alignments of radio 
polarization vectors of the type reported at optical wavelengths by 
Hutsemekers and collaborators, i.e. at cosmological scales. Their 
analysis did not reveal such large-scale alignments at radio wave¬ 
lengths. From this claim, the recognized wavelength dependence of 
the polarization ve ctor alignments has brought the model of axion- 
like particle (e.g . |Das_gt^ llOO ^ IPavez. Cudell & Hutsemekers! 


120081, IXgarwaL Kamal & JainI 201 1[) as the favourite candidate to 

explain alignments at optical wavelengths. This model has however 
been observationally ruled out since it pred icts non-negligible 
circul ar polarization which is not detected Mutsemekers et aU 
I 2 OIOI and lPavez, Cudell & Hutsemekersll^Ol ih . 


Beside this analysis, iTiwari & JainI (1201 3h tested the unifor¬ 
mity of the polarization PA’s considering roughly the same sample. 
They found significant evidences for alignments at distance scale 
of the order of 150 MpcQ. As the correlations are found at differ- 
ent distance scale s, their study d oes not contradic t the analysis of 
Ijoshi et al] ( l2007h . More recently. Tshurtleflj ( l2014li studied the cor¬ 
relation of the PA’s for sources grouped in circular regions of 24° 
radius. While not very significant, he reported indication of PA’s 
alignments in two regions of the sky. 

Despite these analyses which involve different statistical tests 
and different samples which correspond to different cuts of the 
original dataset, the status of polarization PA correlation at radio 
wavelengths in not clear. Moreover, an analysis taking the redshift 
of the sources into account in still missing. The redshift depen¬ 
dence being an important characteristic of the alignments of quasar 
polarization vectors at optical wavelengths, it seemed important to 
us to take it into account in the analysis of the radio sample. Es¬ 
pecially if one seeks the same signature at radio wavelengths as at 
optical wavelengths. Therefore, originally motivated by the redshift 
dependent analysis of this radio sample on the one hand, and by the 
recent and independent confirmation bv IPelgrims & Cudelll l l2014h 
of the polarization alignments at optical wavelengths, on the other 
hand, we devote this study to a careful analysis of the uniformity 
of the polarization PA’s of FSRS’s belonging to the JVAS/CLASS 
8.4-GHz surveys. 

We introduce the data samples which are studied throughout 
this work in Section In Section 13.11 taking the redshift of the 
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^ Attention has to be paid regarding this scale as these authors defined the 
comoving distances assuming a redshift of one for all objects. 
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sources into account, we investigate the polarization PA distribu¬ 
tions of the FSRS’s located in regions of the sky where the op¬ 
tical polarization alignments are the most significant. Stimulated 
by the detection of alignment in one of these regions, we perform 
a complete analysis of the entire dataset in Section [Tl2\ with and 
without accounting for the redshift. Having highlighted significant 
alignment signatures in the sample of quasars (hereafter, QSO), we 
search for their characterization in Section|4l We finally summarize 
our results in Section[3 bring arguments against and for the hypoth¬ 
esis of biases in the dataset and discuss a possible interpretation of 
the data. We conclude in Section that the dataset of the polar¬ 
ization angle measurements of the JVAS/CLASS 8.4-GHz surveys 
are either not exploitable due to instrumental biases or, more excit¬ 
ingly, that we may have pinpointed large-scale alignments at radio 
wavelengths. As we use little known statistical tests throughout this 
work, we found necessary to discuss them. Thereby, the descrip¬ 
tions are brought together in the Appendix in order to keep a fluent 
reading. 


2 DATA SAMPLE 


Table 1. The different samples. 


z 

Object Type 

n 

acronym 


All 

3858 

— 


QSOs 

1450 

QSO 

no 

Radio Sources 

1379 

RS 


Galaxies 

381 

G 


Other Objects 

648 

VO 


All 

1531 

All{z) 


QSOs 

1325 

QSO{z) 

yes 

Radio Sources 

11 

- 


Galaxies 

184 

- 


Other Objects 

11 

- 


Number (n) of the different source species for the objects retrieved in the 
NED database among the sample of 4155 sources with reliable polaiization 
PA measurements, with and without redshift information, z. The last col¬ 
umn contains the acronyms used for the samples analysed in this work. The 
category named Other Objects contains various species with small mem¬ 
bership. 


The JVAS/CLASS 8.4-GHz catalogue is made of the JVAS 
(Jodrell-VLA Astrometric Survey) and the CLASS (Cosmic Lens 
All-Sk y Survey) surveys that were gathered by Ijackson et al.l 
( l2007h to build the largest catalogue of polarization measurements 
of compact radio sources, at that time, pa ying attention t o bias es on 
polarization measurements. We refer to Ijackson et al.l l l2007h and 
references therein for a complete description of the catalogue and 
the surveys. In this catalogue, the total number of object having 
polarization measurements is 12.743 (see the on- line c atalogu a. 
Adopt ing the prescription of Ijackson et al.l (l2007h and Ijoshi et al.l 
(I2007I) . we retain the sources for which the polarized flux is higher 
or equal to 1 mjy in order to keep significant polarization detections 
only and to obtain an unbiased sample. When there is more than 
one object in a radius of 1 arcsec on the sky, we select the object 
with the highest polarized flux. This selection, which also elimi¬ 
nates multiple measurements, leaves us with a sample size of 4265 
obje ctfl Led by a sele ction criterion used at optical wavelengths 
(see iHutsemekerslI 1 9981 for a discussion), we further constrain the 
sample asking that gq ^ 14°, where erg is the error on the polariza¬ 
tion PA. Out of the 4265 sources, 4155 satisfied the criterion. This 
sample, which we call All in the reminder, constitutes the largest 
sample for which we have robust polarization PA measurements 
from the JVAS/CLASS 8.4-GHz surveys. 

Using the NASA Extragalactic Databas^fl (NED), we iden¬ 
tified a total of 3858 sources. We first used the automated mode 
”Near-Object/Position List” with a search radius of 0.1 arcsec. Af¬ 
ter manual clarification of multiple identifications, 3446 objects 
were found. Eor the 709 objects left, we used a search radius of 0.5 
arcsec and found 412 additional sources, after having again manu¬ 
ally clarified the multiple identifications. We stopped the procedure 
at this value of the search radius in order to ensure proper identifi¬ 
cations. Out of the 3858 retrieved objects, 1531 have spectroscopic, 
and thus reliable, measurements of redshift, z. 

The use of NED also leads to the classification of the sources. 


^ http://vizier.u-strasbg.fr/viz-bin/VizieR-3?-source=J/MNRAS/376/371 
® If we only remov e multiple measure ments, we recover the source number 
of 4290 studied by Ijoshi et alj ilOOTl) . We nevertheless choose to add the 
above constraint for an efficient source separation. 
http://ned.ipac.caltech.edu/ 


Table [T] reports the number of sources identified for each class of 
FSRS as well as the number of these sources for which we have 
redshift information. As it can be clearly seen, the QSO’s represent 
86% of the sample with redshift measurements. Hence, analyses 
and results involving samples with redshift information will mainly 
concern those objects. 


3 UNIFORMITY OF POLARIZATION PA’S IN 
JVAS/CLASS 8.4-GHZ SURVEYS 


3.1 Regions of optical polarization alignments 


In 


Hutsraieke^ 1998h. Hutsemekers et al.1 ( l2005h and recently in 


[Pelgrims & Cudelll ( l2014h . specific regions of the sky have been 
highlighted for which polarization PA’s of quasars are found to be 
aligned at optical wavelengths. In the first studies, the two most 
significant regions were identified by eye and were called AI and 
A3. In the latter and independent identification. [Pelgrims & Cudelll 
( I2OI4I) used a more unbiased method and highlighted regions N2 
and S2. While less extended, the latter two regions were consis¬ 
tently found at similar locations in the 3-dimensional space. Here, 
as the sky coverage of the radio surveys and the optical catalogue 
are different, we choose to consider the most extended regions (A1 
and A3) to ensure an overlap as big as possible. Furthermore, those 
regions have been the subject of various studies in the past. They 
are delimited in right ascension, declination and redshift by: 

• Al: 168° < a < 218° ; 5 ^ 50° and 1.0 < ^ < 2.3 

• A3: 320° < a < 360° ; 5 ^ 50° and 0.5 < ^ < 1.5 


where a and 5 refer to the right ascension and the declination of 
the s ources, respective ly. 

Ijoshi et ^ (I 2 OO 7 I 1 addressed the question of uniformity of the 
polarization PA’s of FSRS’s from the JVAS/CLASS 8.4-GHz sur¬ 
veys in these regions and reported no obvious alignment. However, 
they did not introduce the cuts in redshift and thus, only considered 
the windows toward the Al and A3 regions, defined by cuts in right 
ascension and declination only. As the redshift of the sources is an 
important characteristic of the alignments of optical polarization 
PA’s, we perform a new analysis of these regions. It is nevertheless 
important to realize that the sky coverages of the JVAS/CLASS 


(c) 2015 RAS, MNRAS POQ.fTlflTI 













































Polarization Alignments of Radio Quasars 3 


FSRS’s and the studied quasars at optical wavelengths are differ¬ 
ent. In particular, the surveys at radio wavelengths do not contain 
data at (5 < 0°. 


To study the polarization PA distributio ns of the A1 and 
A3 re gion s, we use the Haw ley-Peebles test dHawlev & Peeble^ 
I 1975 I and Godlowski ^12h and the density test introduced in 
dPelgrims & Cudellll2014ll . We refer to Sections lAll and IA2I for a 
description of these statistical tests. Considering samples of sources 
that are in a small region of the sky, these tests return the probabil¬ 
ity that the observed distribution of PA’s is random and define the 
mean PA’s (9 and 9pc, resp.) which are relevant only in the case of 
non-uniformity. 

The Hawley-Peebles test analyses PA histograms. The num¬ 
ber of bins is a free parameter. We decide to use 18 bins of 10° 
each, spanning the range 0° — 180°. This choice does not maxi¬ 
mize the reported probabilities but is somehow justified by the fact 
that a bin width of 10° corresponds approximately to twice of the 
mean of the errors of the polarization PA’s under study. 

The density test involves spherical caps of equal area. This 
area is fixed by the angular aperture of the cap (2rj) which is the 
only parameter of the test. As discussed in the Appendix, it is useful 
to investigate a large range of values for rj. Therefore, we explore 
here the range 2° — 90° for rj and we report the probabilities Pmin 
and corresponding to the value of rj for which is the smallest. 

Results of the tests applied to the sub-samples extracted from 
the All{z) and the QSO{z) samples are shown in Table The 
hypothesis of uniformity of the polarization orientations is rejected 
at the level of at least 95% in the A3 region but not in the A1 re- 
gioifl The discov ery of this alignme nt is very intriguing, especially 
given the claim of ijoshi et alj 1 1200'^ . To reconcile the analyses, we 
extracted the A3 window from the sample of 4155 FSRS’s. Ap¬ 
plying our tests to this sub-sample of 385 sources, we did not find 
any eviden c e for alignment and thus confirm the negative result of 
Ijoshi et aP ( l2007h (see Tablej^and Section 133] for further related 
discussions). 


As a result, we contradict the claim of jjoshi et aP ( l2007t) 
which states that no alignment is present at radio wavelengths (8.4 
GHz) inside the regions where the optical polarization vectors are 
found to be coherently oriented. This contradiction is likely due to 
the redshift cuts. 


It is of interest to find out whether the alignment tendency ob¬ 
served in this region is an isolated structure inside the sample of 
FSRS’s or if it is part of a major trend which was not recognized 
by earlier works. To this end, we shall address the question of uni¬ 
formity of the polarization PA’s for the complete JVAS/CLASS 8.4- 
GHz surveys, without restriction on the sky location, taking redshift 
into account and considering the subdivision of the sample into the 
different species. 


3.2 Full sky coverage 

To study the uniformity of polarization angle distributions for 
samples of sparse and non-uniformly scattered so urces on the 
celesti a l sphere, tools have been d eveloped by Hutsemekep] 
(Il998l) . jjain, Narain & Saralal (|2004|) . iHutsemekers et al.l ( l2005h . 


® At first glance, the fact that we found alignment in the A3 region but not 
in the A1 region could be caused by a bad spatial overlap between radio 
and optical data in the A1 region. But the overlapping is not better in the 
A3 region. 


jjoshi et aP (|2007|) and bv jPelgrims & Cudell (2014). The so-called 
S and Z test s, estab lished by iHutsemekerd d 19981) and modified 
by jjain et all ( |2004|) to obtain coordinate-invariant statistics, are 
appropriate to assess the probability that the distributions of po¬ 
larization PA’s of local groups are due to statistical fluctuations 
considering the overall sample. We use these tests in this section. 
The intrinsically coordina te-invariant statistical test introduced in 
fPelgrims & CudelH )|2014|) being more useful for the characteriza¬ 
tion of correlations is used in Section |4] We do not use the other 
statistical tests as they are coordinate-dependent, this dependence 
growing with the angular distance between sources. They are thus 
not adequate to test the uniformity of the PA distribution over large 
scales. 

The coordinate-invariant S and Z tests are extensively dis¬ 
cussed in Section |A3] These nearest-neighbour tests compute the 
probability that the polarization PA’s are uniformly distributed in 
spatially defined groups of objects, making use of Monte Carlo 
simulations. For each realization, the PA’s are reshuffled among the 
sources of the entire sample and a statistics is computed for each 
group of n„ nearest neighbours. The percentage of Monte Carlo 
simulations having an average statistic {Sd or Zc) as extreme as 
the one of the data defines the significance level (S.L.) of the test, 
i.e., the probability that the observed PA correlations inside groups 
can be attributed to statistical fluctuations in the entire sample. To 
evaluate the S.L., we set the number of random simulations to 1000, 
except contraindication. For the samples of Table [T] we explore a 
wide range of values of the parameter riv (see Section lA3.4l for the 
motivation of doing this). We span the range 4 — 400 (with steps 
of 2 and 20 for ranges 4 — 18 and 20 — 400, respectively), except 
for the sample of galaxies where we stop at n„ = 200 for obvious 


We first consider samples for which reliable redshift measure¬ 
ments are available. For these samples, both 2- and 3-dimensional 
analyses are applied, i.e. defining nearest-neighbour groups on the 
celestial sphere or in the 3D comoving space respectively. We then 
turn to the 2-dimensional analysis of the samples of Table [T] which 
are not constrained by redshift. For convenience, we give in Ta¬ 
ble |4] a summary of the results of these two tests applied to all the 
samples. 


3.2.1 Samples with redshift measurements 

As already mentioned, the sample of 1531 sources for which re¬ 
liable redshift measurement are available is composed at 86% by 
QSO’s. The second important population of this sample is that of 
galaxies. The redshift distributions of these samples are shown in 
Fig.m Of course, the redshift distributions of the sample QSO{z) 
and the one of galaxies do not follow the same trend. 

3.2.1.1 3-dimensional analysis. We discuss our results in term 
of a typical comoving distance £ instead of the free parameter riv 
of the statistical tests. This typical scale is defined as being the me¬ 
dian of the comoving distances between each central object and its 
n^-th nearest neighbour, the median being evaluated over the full 
sample under consideration. The line-of-sight comoving distances 
are computed assuming a flat Universe with the cosmological pa¬ 
rameters: Gm = 0.31 and Hq = 68 kms“^Mpc“^. We show in 
Fig. ID (Top) the relation between the parameter value n„ and the 
typical comoving distance for the samples we analyse in 3 dimen¬ 
sions. We applied the statistical tests to the sample of 1531 objects 
and to the subcategory of QSO, namely All{z) and QSO{z). We 
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Table 2. Test results in the A1 and A3 regions. 



A1 

region 

A3 region 


All(z) 

QSO{z) 

AU(z) 

QSO(z) 

n 

141 

139 

50 

45 

Fhp (%) 

95.43 

93.15 

4.19 

0.96 

e(°) 

- 

- 

64 

62 

Pmin (^) 

^ 4 

^ 5 

0.01 

7.2 lO-'^ 

v{°) 

20 - 90 

20 - 90 

58 

56 

p'^ (%) [risim] 

^ 13 [102] 

^ 10 [102] 

0.07 [510-*] 

4.010-3 [510-*] 

dpc (°) 

- 

- 

68 

67 


Results of the Hawley-Peebles and the density tests performed on the sub-samples corre¬ 
sponding to the A1 and A3 regions of optical polarization alignments. Sub-samples are ob¬ 
tained from both the All(z) and the QSO{z) samples, n is the size of the sub-samples, Php 
is the probability given by the Hawley-Peebles test that the PA’s are drawn from a uniform 
parent distribution and 0 is the mean polarization PA returned by this method, pmin is the 
local probability obtained with the density test for the half-aperture angle rj to which corre¬ 
sponds the minimum global probability computed with risim random simulations. The 
mean angle Opc is computed as explained in the text. Probabilities are given in percent. 



z — redsliift 


Figure 1. Redshift distributions of the sample All(z) and its sub-samples 
of QSO’s and galaxies. 


also considered the high redshift part of the latter, imposing z > 1. 
This restricted suh-sample is populated by 894 sources and is de¬ 
noted QSO{z > 1). 

Results of the S and Z tests are shown in Fig. (Middle and 
Bottom, respectively). We did not find any significant evidence 
(S.L. < 5%) over a wide range of value of n„ (or £ ) for align¬ 
ment of the polarization PA’s in the sample All{z) and QSO{z). 
However a redshift dependence is possibly detected with the Z test 
as suggested in Fig. (Bottom). Indeed, in the high redshift QSO 
sample, correlations of polarization PA’s of sources inside groups 
of typical comoving radius of about 2 Gpc show a probability less 
than 1% of being due to statistical fluctuations. 

3.2.1.2 2-dimensional analysis. For the 2-dimensional analysis, 
the radial coordinate of the sources is fixed to r = 1, even though 
redshift measurements are available. We discuss our results in term 
of the typical angular separation The latter is defined as being the 
median of the angular separation between each object of the sample 
and its nearest neighbour. Fig. (Top) shows the relation 

between ^ and n„ for the three samples. 

We show the dependence of the S.L. on the typical angular 


separation ^ in Fig.|^ Significant correlations (S.L. < 5% over a 
wide range of ^ value) of the polarization PA’s inside groups is ob¬ 
served for the three samples (AU{z), QSO{z) and QSO{z > 1)) 
although the minima occur at different typical angular separations. 
Indeed, for the S test, the sample All{z) shows its minimum S.L. 
at 0.3 % for ^ « 23° and QSO{z) shows a small dip for the range 
of ^ « 8 — 26° with a minimum S.L. = 1.2% for ^ = 18°. 
The high redshift part of the QSO sub-sample {QSO{z > 1)) ex¬ 
hibits values of the S.L. below 1% for smaller angular separation 
(^ < 10°). These features are somehow confirmed by the Z test as 
seen from Fig. [3 (Bottom). For this test, the minimum S.L. value 
of the sample QSO{z) is found to be as low as 0.3 % for ^ « 23° 
and the sample QSO{z > 1) shows S.L. below 1% for ^ ^ 10° 
with an additional dip around ^ = 34°. 

3.2.2 Full samples with different object types 

Analysing samples with redshift measurements in two dimensions, 
we have found significant correlations (S.L. < 5% with minima 
< 1%). It is therefore interesting to also perform the 2-dimensional 
analysis on the other samples of Table [T] i.e. on the samples not 
restricted by the availability of the redshift of the sources. We thus 
consider the sample of 4155 sources as well as its four sub-samples 
with different object types. We show in Fig. |4] (Top) the relations 
between n„ and ^ for these five samples. The results of the S and Z 
tests are shown in Fig.|4](Middle and Bottom respectively). 

For small values of n„ (from 6 to 10), we found indications 
of alignments in the sample All as the S and Z tests return S.L. 
values at the percent level (1.1% and 1.2%). These indications 
of alignments are re miniscent of the correlations highlighted by 
iTiwari & JainI ( l2013h at the scale of ~ 150 Mpc. The reasons why 
we found less convincing correlations are likely that (/) we consider 
a different sampl^fl and (ii) that we use a different statistics. 

For large values of n„, alignments are detected with S.L. be¬ 
low 5% over a wide range of ^ only for the sample QSO. The 

® When we built the dataset fr om the JVAS/CLASS catalogue, we removed 
duplicate measurements while iTiwari & JainI J2013h did not (Jain (private 
communication, 2015)). 
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Figure 2. 3D analysis of the samples with redshift measurements. Top: Re¬ 
lation between the parameter riv and the typical comoving separation C in 
Mpc for the samples ylZZ( 2 ), < 350 ( 2 :) and ( 350 ( 2 : > 1). Middle and Bot¬ 
tom: Significance level obtained with the S and Z tests (resp.) as a function 
of the typical comoving distance C for the three samples. The 5% and 1% 
S.L. are indicated. 


Figure 3. 2D analysis of the samples with redshift measurements. Top: Re¬ 
lation between the parameter rt„ and the typical angular separation ^ in 
degree for the samples All{z), QSO{z) and QSO{z > 1). Middle and 
Bottom: Significance level obtained with the S and Z tests (resp.) as a func¬ 
tion of the typical angular separation ^ for the three samples. The 5% and 
1% S.L. are indicated. 


S.L. of the S test applied to QSO exhibits a dip for n„ between 
40 and 140, reaching the value of 0.7 % for = 60 and 80. 
The range of typical angular separations involved in these cor¬ 
relations is ^ « 12° — 24°, with stronger correlations for the 
range 14.5° — 17.5°. The Z test exhibits a large dip for the range 
n„ = 40 — 200 with the minimum S.L. of 0.12 % for n„ = 140, 
implying correlations at ^ « 24°. Those correlations of polariza¬ 
tion PA’s involve QSO’s separated by large distances on the celes¬ 


tial sphere and confirm the detection made in the sample QSO{z) 
(with redshift measurements) in Section [3.2. II 

For such angular scales, the distributions of the polarization 
PA’s of the other samples {All, RS, G and VO) are in good agree¬ 
ment with the hypothesis of uniformity. Let us emphasize that these 
large-scale correlations are not observed for the category of radio 
sources (RS), even though the sample size is comparable to the one 
of QSO (see Table[T]l. 
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Figure 4. 2D analysis of the full samples. Top: Relation between the pa¬ 
rameter riv and the typical angular separation ^ for the samples All, QSO, 
RS, G and yO of Table[T] Middle and Bottom: Significance level obtained 
with the S and Z tests (resp.) as a function of the typical angular separation 
^ for the five samples. Note that the S.L. value of the sample QSO for 
^ Si 24° (n„ = 140) has been computed with 10^ random simulations. 


(l2007h . As already mentioned, when these authors analysed the so- 
called A1 and A3 regions, they did not constrain their sample with 
regard to the redshift, which is an important characteristic of the 
optical polarization alignments. They applied their Nearest Neigh¬ 
bour Test to the full sample restricted to the sky window of the 
A3 region, i.e. introducing cuts in right ascensio n and declination 
only (see Section 3.3 and 5 of Ijoshi et alj|200% . The first possi¬ 
ble cause of divergent results is a redshift dependence of the align¬ 
ments at radio wavelengths as at optical wavelengths, so that taking 
all the sources in the window regardless of their redshift blurs the 
alignment. Alternatively, the comparison of the last two columns 
of Table 12 suggests that alignments are more pronounced for QSO 
compared to the other types of objects. 

In order to test these two scenarios, we performed an analysis 
of the samples obtained by imposing the A3 window cut on the dif¬ 
ferent samples of Table [T] Results of the Hawley-Peebles test and 
the density test of Pelgrims and Cudell are shown in Table [2 Cor¬ 
relations between polarization PA’s are observed when we consider 
the A3 window cut of the samples QSO, QSO{z) and All{z), but 
no deviation from uniformity is detecte d for the A3 wind ow cut 
of All, in agreement with the result of Ijoshi et al.l (l2007h . Com¬ 
parison of the last two columns of Table [2 teaches us that adding 
the other species to the QSO sample completely blurs the align¬ 
ments. This simple observation argues for the scenario in which the 
species selection is at the origin of the detection of the correlations. 
This scenario is reinforced when we consider the A3 window cut 
of the sample RS. For this sub-sample of 138 objects, we found 
that the distribution of the polarization PA’s is in agreement with 
the hypothesis of uniformity. Comparison of Tables |2 and [2 does 
not allow us to conclude on a possible redshift dependence of the 
polarization alignments. This is partially due to the lack of redshift 
information for species other than QSO’s. 

In the full sample we highlighted alignments involving 
sources separated by typical angular scales of about 20°. For the 
samples with redshift measurements, the large-scale correlations 
are observed to be more significant with the 2-dimensional analy¬ 
sis than with the 3-dimensional one. Considering samples that are 
not limited by the redshift availability, we also pinpointed that the 
large-scale correlations mainly concern the category of QSO as it 
was already suggested during the study of the A3 window in Sec¬ 
tion |2l] 

As a conclusion, we find that the polarization PA’s of the 
JVAS/CLASS 8.4-GFIz surveys show correlations in groups of 
QSO’s with an angular radius of about 20°. The significance level 
at which these correlations can be attributed to statistical fluctua¬ 
tions in the sample of QSO is found to be as low as ~ 0.1% for 
^ « 24° with the Z test. 


3.3 Preliminary results 

So far, studying the polarization PA’s of different samples drawn 
from the JVAS/CLASS 8.4-GHz surveys, we have found significant 
alignments in some of these samples; first, in one of the regions 
where the optical polarization vectors were found to be aligned and 
second, in the QSO all-sky survey. 

Regarding the A3 region, a few reasons might lea d to the dif- 
ferences between our findings and the conclusions of Ijoshi et al.l 


4 IDENTIFICATION OF REGIONS OF ALIGNED 
POLARIZATIONS 

For the correlations highlighted in the previous section, it would be 
of interest to figure out if the alignments detected at typical scales 
of ^ « 15° — 25° are due to a global trend across the whole sky 
coverage of the survey or if they are prominent in some regions of 
the sky, as it seems to be the case at optical wavelengths. 

To this end, we proceed to the identification of the groups of 
sources with distributions of polarization PA’s that show significant 
departure from uniformity. The fact that these groups are clustered 
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Table 3. Test results with the A3 window cut on QSO{z), All(z), QSO and 
All samples. 



QSO(z) 

A3 window cut 

All{z) 

on 

QSO 

All 

n 

100 

115 

114 

385 

Fhp (%) 

0.36 

1.27 

1.19 

29.6 

en 

68 

72 

68 

- 

Pmin (%) 

7.0 10-“* 

1.910-3 

5.910-3 

0.14 

h(°) 

52 

52 

52 

68 

(%) [riaim] 

0.02 [510'*] 

0.01 [10^] 

5.410-2 [510^1 

1.5 [103] 

®PC (°) 

69 

70 

68 

59 


Same as Table|2]but for the sub-samples obtained by application of the A3 win¬ 
dow on the samples of 1325 QSO’s with redshift, 1531 sources with redshift, 
1450 QSO’s regardless of the redshift information, and 4155 flat-spectrum radio 
sources. 


Table 4. Summary of S and Z test results. 




S 


Z 



3D 

min(S.L.) (%) 

£ (Gpc) 


min(S.L.) (%) 

C (Gpc) 

Uv 

All(z) 

- 

— 

— 

— 

— 

— 

QSO(z) 

- 

- 

- 

- 

- 

- 

QSO(z > 1) 

- 

- 

- 

0.6 

~ 1.7 

40 

2D 

min(S.L.) (%) 

^(°) 

riv 

min(S.L.) (%) 

«(°) 

riv 

All{z) 

0.3 

~ 23 

140 

1.0 

~ 23 

140 

QSO(z) 

1.2 

~ 18 

80 

0.3 

~ 23 

120 

QSO{z > 1) 

0.5 

~ 10 

18 

0.6/0.5 

10/34 

18/160 

All 

1.1 

3-4 

8-10 

1.2 

~ 3 

6-8 

QSO 

0.7 

14.5 - 17.5 

60 - 80 

0 .12* 

~ 24 

140 

RS 

G 

VO 

1.3 

3 

4 

1.1 

~ 3 

4 

- 

- 

- 

3.0 

40 

160 


Summary of the application of the S and Z statistical tests to all samples of Table[T] For each test, we report 
the value of the minimum S.L. with the corresponding n„ parameter and its attached typical scale {C or ^ 
for the 3- or 2-dimensional analysis, resp.). We only show results when the S.L. of the sample is found to 
be below the threshold of 5% for a wide range of Uv All S.L. have been evaluated with 1000 Monte Carlo 
simulations except the smallest one (marked by an asterisk) for which we had to use 10000 simulations. 


in space or not tells us whether the correlations of polarization ori¬ 
entations are due to well localized objects or to a general trend. 
This identification can, inter alia, be achieved with the help of the 
S and Z tests. For clarity, we give the details for the S test. Also, 
note that we limit our search to the 2-dimensional analysis of the 
sample QSO since it revealed the most convincing evidence of de¬ 
parture from uniformity with confidence level higher than 99% for 
the range of f « 14.5° - 17.5°. 

In Section 13.21 local statistics Si were computed for each 
nearest-neighbour group and these statistics have been computed 
for each simulated dataset. We attribute to each central source i the 
quantity Si which tells how much the corresponding group of near¬ 
est neighbours contributes to the global statistics Sd ■ This quantity 
is defined as Si = {{Si} — Sf jloi, where S* is the statistics 
obtained for the observed dataset (see EQS. IAll and lA2t and where 
{Si) and Gi are the mean and the standard deviation of this statistics 
evaluated over the whole set of simulations. The bigger the value 
of Si, the more the group contributes to So- A group of nearest- 


neighbour objects is considered as contributing significantly to Sd 
if Si ^ Sc, for a given threshold Sc. 

Of course we shall search for the identification of the most 
significant groups, i.e. consider the Si quantities computed with the 
parameter riv chosen such that Sd is the smallest (see Table|4l(. To 
visualize the sky location of the most significant groups we pro¬ 
duce maps which highlight their central sources. Note that these 
maps do not critically depend on t he choice of n„. Th ese maps are 
equal area Schmidt projection (e.g. lFisher et al.ll987h of the north¬ 
ern hemisphere (in equatorial coordinates). This choice is suitable 
for the considered dataset as it covers only positive declinations. 
We also plot (the grey bold lines) the limits of the A1 and A3 win¬ 
dows defined in Section [TT] Let us insist on the fact that only the 
northern part of these limits are shown. The A1 and A3 windows 
defined from the analysis at optical wavelengths extend to the South 
hemisphere which is not displayed here. 

The identification map corresponding to the S test (in 2D) for 
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Figure 5. Identification map for the sample QSO using the S test in 2D. 
Parameters are fixed to = 80 (see Fig. |4] (Middle)) and two thresh¬ 
old values: Sc = 2.5 (lighter points) and Sc = 3.0 (darker points), re¬ 
spectively cyan and light blue (color on-line). Identification maps are equal 
area Schmidt projection of equatorial coordinates. Only the equatorial north 
hemisphere is displayed with the north pole at the centre of the map. Grey 
circles are parallels of declinations 0°, 30° and 60° and grey diagonals are 
meridians of right ascensions being multiple of 45°. The curved dashed line 
is the Galactic equator, the North and the South Galactic caps being respec¬ 
tively above and below the line. Grey bold lines are northern boundaries of 
the A1 and A3 regions of optical polarization alignments (see text). Small 
black dots are locations of the 1450 sources of the sample. Highlighted 
sources are objects for which coiTesponding neighbours show a polariza¬ 
tion PA distribution that is unlikely due to chance. 


the sample QSO with the parameters = 80 and Sc = 2.5 is 
shown in Fig. |5] As one can see, the highlighted central sources 
cluster in three or four groups along with other more sparse and/or 
isolated locations. When pushing Sc up to 3 (darker points), only 
the cluster with right ascension a ~ 206° and declination 5 ~ 38° 
remains. Following this analysis, it is likely that the significant de¬ 
parture from uniformity in this sample is due to polarization align¬ 
ments in a few groups of QSO’s. It is intriguing that two of them 
are found in the A1 and A3 windows. 

In order to put the latter identification of aligned regions on 
stronger grounds, we may use other tests. The Z test also reveals 
significant non-uniformity. For the QSO sample and the parameter 
value = 140, it leads to the map shown in Fig. |6] which is 
in relatively good agreement with Fig. although it shows more 
scattered clusters. 

Also, we find relevant to p roceed to a complementar y identi¬ 
fication using the density test of lPelgrims & Cuddl l l2014h . As this 
test does not depend explicitly on the number of nearest neighbours 
(while it is encapsulated within the statistics), local groups can be 
defined by a physical angular scale, denoted fl. In order to carry 
out an identification as close as possible to those produced with the 
S and the Z tests, we found necessary to split each sample in its two 
Galactic hemispheres to determine the physical scale at which local 
groups have to be defined. Indeed, the density of the data points in 
the North Galactic cap and the South Galactic cap are different. In 
the sample QSO, the typical angular separation corresponding to 
= 80 is 5 ~ 17° and corresponding to = 140 is ^ « 24° 
(cf. Fig.|4](Top)). However, by splitting the sample in its northern 


Figure 6. Identification map for the sample QSO using the Z test in 2D. 
Parameters are fixed to = 140 (see Fig. |4] (Bottom)) and Sc = 1.65 
two threshold values: Sc = 1.65 (lighter points) and Sc = 1.75 (darker 
points), respectively cyan and light blue (color on-line). Please note that the 
thresholds Sc for the S and Z tests do not refer to the same quantities and 
have thus different values. 


and southern Galactic parts, we obtain ^jv « 16° and ~ 23° for 
riv = 80 and ^jv « 21° and ~ 33° for Uv = 140, respectively. 
Given these values, we decided to define local groups in 2 dimen¬ 
sions with angular scales Hjv = 20° for the North and Qs = 30° 
for the South. 

As we search for the characterization of the polarization PA 
distribution of each group taken as a whole, we shall not investi¬ 
gate values of rj (the free parameter of the method) below the re¬ 
spective angular separation, i.e. below the imposed angular scales. 
We arbitrarily chose rj — 40° and rj = 50° for the North and the 
South, respectively. The identification map computed with these 
parameters is shown in Fig.|7] We checked the robustness of the 
map with other pairs of values such as (Hat, fls) = (15°, 25°) 
and (25°, 35°). We also checked the stability of our results using 
other values of t). Note that we did not search for the optimal value 
of r), i.e. the one which would give the lowest probabilities, as the 
method undergoes edge effects. We rather spanned the range of 20° 
to 60° with step of 5° and found consistent maps. 

Although a close examination shows discrepancies in the pre¬ 
cise locations of central sources of neighbouring groups, compar¬ 
ison of maps presented in Figs. [3|6] and |3 reveals relatively good 
agreement, especially for the cluster at (a, 5) ~ (206°, 38°). 

In order to define more precisely the limits of regions of po¬ 
larization alignment, we proceed as follows. To each central source 
corresponds a group of nearest objects (defined via the parame¬ 
ter Uv or H). A highlighted central source is said to form a cluster 
along with (an)other highlighted source(s) if it belongs to the group 
of nearest objects of the latter. A central source is discarded from 
a cluster if it is not in the neighbourhood of a sufficient percentage 
of central sources forming this cluster (e.g. ~ 60%). Reproducing 
this test for all highlighted objects, we end up with identification 
of independent clusters. We finally add to the cluster the nearest 
neighbouring objects of each central sources, paying attention to 
duplication. Although this procedure is rudimentary, it is sufficient 
for our goal. 
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0 


Figure 7. Identification map for the sample QSO with the density test. 
Parameters are fixed as = 20 ° and rj = 40 ° for the North Galactic 
cap and Q = 30 ° and r} = 50 ° for the South Galactic cap. Highlighted 
sources correspond to groups showing Pmin ^ 10“® (lighter points) and 
Pmin ^ 3 10 “^ (darker points), respectively cyan and light blue (color 
on-line) 


Table 5. Identified regions from the 2-dimensional analysis of the sample 
QSO with the S, Z and density tests. 



n 

(a, <5)cm (°) 

«(°) 

^•max ( ) 

Php (%) 

en 

RNl 

108 

(163, 12) 

12 

21 

0.45 

131 

RN2 

191 

(206, 38) 

14 

25 

1.17 

42 

RSI 

116 

(340, 18) 

15 

25.2 

1.45 

57 


Regions are intersections of those given by each test (see text). The two first 
lines are for the regions located in the North Galactic cap and the third is 
for the region of the South part. They are named RN1, RN2 and RS1, re¬ 
spectively. n is the number of members belonging to the region, [a, (5)cm 
refers to the position of the normalized vectorial sum of the source loca¬ 
tions of the region, ^ and ^max are the mean and the maximum value of the 
angular separations of sources to (a, <5)cm- S are the results of 

the Hawley-Peebles test. 

We thus end up with three regions for each of the three tests. 
We decide to define our final regions as being the intersection of 
the regions from the different tests. We report the final regions of 
alignments in Table with some characteristics and the result of 
the application of the Hawley-Peebles test on their polarization PA 
distributions. 

As a result, we identified three well-defined regions of the sky 
in which QSO’s show coherently oriented polarization vectors. Two 
of these regions are located in the North Galactic hemisphere of 
the sky and one toward the South. Considering the southern cap, 
it is worth to remark that more than 85% of the sources of the 
sub-sample identified here belong to the A3 window defined from 
the region of optic al alignment discovered by Hutsemekers et al. 
( 119981120^ l2()05h . To the North, the identified regions are some¬ 
how located at the edges of the Al window of optical alignment, 
one outside at low declination and the other inside at high decli¬ 
nation. We call them RN 1 and RN2, respectively. It is again worth 


mentioning that more than 70% of the sources of the RN2 sub¬ 
sample identified here belong to the Al window. 

It is remarkable that our region RN2 coincides with the main 
aligne d cluster resulting from the independent analysis of lShurtl^ 
( I 2 OI 4 I) . Consistent with our previous results, we report a stronger 
alignment than he did as we only consider the species of QSO. 

In order to visualize the alignment patterns, we show in Fig. [8] 
the equatorial-coordinate maps of the normalized polarization vec¬ 
tors of the identified regions along with their corresponding polar¬ 
ization PA histogram. Some structures can be spotted out by eye. 
This is better seen in the region RNl (see Fig. [8] (Top)). The sta¬ 
tistical tests used throughout this study do not allow us to better 
search and characterized such structures. This task is far beyond 
the scope of this paper and would request dedicated algorithms to 
compute the likelihood of structures of aligned polarization vectors 
in a random sample. 


5 INTERPRETING THE RESULTS 

So far we found that the polarization vectors of QSO’s being in 
groups having angular radius of about 20° have correlated orien¬ 
tations. We showed that these groups cluster in three independent 
regions of the sky and that to each of these corresponds a different 
preferred polarization PA. 

5.1 Are the data contaminated by systematics? 

The preferred angles for the two northern regions are found to have 
values close to 45° and 135°. These values, being very particular, 
lead us to consider the possibility that the correlations we found 
are d ue to biases in the dataset (see iBattve. Browne & JacksonI 
I 2 OO 8 I) . This hypothesis, being a priori difficult to reconcile with 
the local character of the alignment features, could potentially ex¬ 
plain that they are better detected with the 2-dimensional analy¬ 
sis t han with the 3-dimen sion al one. In th i s sens e, and contrary to 
what I Jackson et all ( l2007h and ijoshi et al.l ( |2007|) claimed, we also 
find evidence for a global non-uniformity inside the polarization 
dataset. Using the Hawley-Peebles test, the probability that the dis¬ 
tribution of the 4155 objects is uniform is found to be Php = 2.7% 
(with 9 ~ 51°). This non-uniformity of the overall polarization 
distribution of the sample All argues for the hypothesis of a biased 
dataset. However, consistently with our previous results, this non¬ 
uniformity is found to come from the sub-category of QSO as we 
find Php = 1.1% (with 9 ~ 57°) for this sample and that remov¬ 
ing the QSO’s from the sample All leads to Php = 44.59CI This 
result together with the previous evidence for alignment of QSO’s 
and not for the other species is awkward to reconcile with an ob¬ 
servational bias, as there is no reason for a contamination of the 
polarization data for the species of QSO and not for the others, as 
we shall see. 

Comparing the properties of the samples QSO and RS 
(which have a comparable number of objects), we note some differ¬ 
ences. As illustrated by the Fig.|3 their polarization characteristics 
at radio wavelengths do not follow the same parent distribution. 
However, we do not find obvious reasons why the sample QSO 
would be more affected by observational biases than the sample 
RS as the QSO sample shows higher total and polarized flux. 

^ The polai'ization PA distribution of the sample RS is also in good agree¬ 
ment with uniformity (Php = 79.1%). 
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Figure 8. Maps of polarization vectors in the identified regions along with 
their corresponding polarization PA histograms. Polarization vectors are 
normalized to the same length in each map. Top and Middle: the two re¬ 
gions in the North Galactic hemisphere at low and high (equatorial) dec¬ 
lination, respectively. Bottom: the region identified in fhe South Galactic 
hemisphere. Properties of these regions are given in Table|^ 



Degree of Linear Polarization (%) 


Figure 9. Polarized flux (Top) and degree of linear polarization (Bottom) 
of the samples QSO and RS. A two-sample Kolmogorov-Smimov test 
reveals that the polarized flux, as well as the degree of linear polarization, 
of the two samples have a probability much below 1% to be drawn from the 
same underlying parent distribution. 


The distribution on the sky of the two samples is also different. 
While the sample QSO is almost homogeneously distributed over 
the sky, the sample RS is far from being so. In order to test the pos¬ 
sibility that it is the difference of sky distribution that is responsible 
for the detection of alignments for QSO and not for RS, we test 
the uniformity of the polarization PA of the RS’s belonging to the 
regions of alignments of Table The overlap is very poor for the 
regions in the North Galactic hemisphere: only 32 RS’s are found 
in each of the RN1 and RN2 regions. To the South, however, there 
are 165 RS’s in the RSI region. The Hawley-Peebles test does not 
reveal departure from uniformity, neither taking RS’s alone or mix¬ 
ing them with the QSO’s of this regiotjf]. Therefore, while the bad 
overlap between the QSO and RS samples in the northern regions 
could explain the difference in the alignment detection, this is not 
the case to the South. We thus conclude that the difference of sky 
distribution is unlikely responsible of this difference. 

Similarly to an instrumental bias, contamination by fore¬ 
ground polarization would affect more strongly the sample of RS’s 
than the one of QSO’s as the polarized flux is globally smaller for 
RS (see Fig. (Top)). This is again in contradiction with what is 
observed. The contamination by foreground polarization is thus un¬ 
likely responsible for the observed correlations of the polarization 
PA’S of QSO’s. 


5.2 Are the polarization alignments real? 


As polarization is usually 
of the object (e. 


correlated to 
Saikia & Salted 


the 


morpho l ogical 
Listed 200 iL 


axis ot the obiect (e.g., |S aikia_&.SjiltM LmeiHzU^ 

IPollack, Taylor & Zavalar2003LISmith et al.l2004l and !Mariij2014h 
there might be real differences between the classes of QSO and 
RS. Indeed, the core dominated FSRS’s are predominantly quasars 


This was expected from the analysis of the A3 window (see Section iTSl . 
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or BL Lac objects in which the jet is oriented close to the line- 
of-sight (see I Jackson et al.ll2007h . The majority of the sources be¬ 
longing to the sample RS is thus expected to be BL Lac objects 
which are thought to be viewed at very small angles to the line- 
of-sight. Consequently, they are expected to show rapid variations 
(< 2 years) of their polarization PA and to be more strongly po¬ 
larized than quasars, which is what we observe in Fig. [9] (Bot¬ 
tom). Also, for this class of object, no net correlation between 
the jet orientation a nd the polarization PA has been reported (e.g. 
jPollack et alj 120031 and references therein). These observational 
facts could bring an explanation to the absence of alignment signa¬ 
tures for the sample of RS’s within the hypothesis that polarization 
alignments reflect morphological-axis alignments of the sources, 
as supported by the recent discovery at optical wavelengths of 
such correl ation for one of the mos t largest known quasar group 
at 2 1.3 jHutsemekers et alll2014h . The latter hypothesis is also 

reinforced by the discovery of large-scale alignments of the jet po- 
sition angles of active galac tic nuclei in the ELAIS N1 field by 
Ijagannathan & Tavlotl ( 12014) and ijagannathani ( l2014h . The fact that 
radio and optical alignments are found in the same parts of the sky 
also supports a real effect. In this framework, one would have to 
compare the alignment patterns observed at optical wavelengths 
with these at radio wavelengths. However, given the bad overlap 
between the sky coverage of the radio and optical catalogues, a de¬ 
tailed comparison is not straightforward. 


6 CONCLUSION 

We tested the hypothesis that the polarization position angles 
are randomly distributed among the FSR S’s contained in th e 
JVAS/CLASS 8.4-GHz surveys presented in ijackson et al.l l l2007h . 
We performed the analysis in two and three dimensions, accounting 
for the distribution of the sources on the sky (2D) and additionally 
for their line-of-sight comoving distances (3D). The polarization 
orientations of radio quasars (QSO sample) show low probabili¬ 
ties to be consistent with the hypothesis of randomness. This de¬ 
parture from uniformity is likely due to correlations of polariza¬ 
tion vectors of QSO’s in groups of angular radius of about 20°. 
A basic identification procedure showed that these groups clus¬ 
ter in three distinct regions of the sky. Two of them fall in the 
A1 and A3 windows of the sky wher e optical polarizati on align¬ 
ments were found (Hutsemekers et al. 1 1 9981 . [200 iL I 2 OO 5 I) . Among 
sources in the JVAS/CLASS sample, only the class of QSO ex¬ 
hibits such large-scale correlations. If real, such alignments at radio 
wavelengths would support the interpret ation of alignments at opti - 
cal wavelengths by spin-axis alignments dHutsemekers et al.l2014h . 
However, our findings prove to be difficult to interpret either as re¬ 
sulting of biases in the dataset or as being the signature of a phys¬ 
ical effect. Indeed, one can find arguments for and against each 
scenario. Among them, the fact that the alignments are more pro¬ 
nounced in 2D than in 3D and that the mean PA’s are multiple of 
45° in some regions would suggest a biased dataset whereas the 
detection of alignments for one class of object but not for the oth¬ 
ers and the clustering of aligned sources in a few regions of the sky 
consistent with those found at optical wavelengths might be seen 
as the signature of a physical effect. 

In conclusion, we highlighted correlations between the radio¬ 
polarization vectors of quasars which could demonstrate the pres¬ 
ence of the same kind of alignment effect as seen at optical wave¬ 
lengths, or alternatively, which could demonstrate that the ra¬ 
dio polarization catalogue is affected by observational biases and 


thus cannot be used to study the polarization o rientations of flat - 
spectrum radio sources. Therefore, the claim bv lloshi et al.l l l2007h 
stating that, at radio wavelengths, there is no alignment signature 
of polarization vectors on cosmological scales of the type found at 
optical wavelengths should be seen with caution. 

More data are clearly needed to assess either the reality of po¬ 
larization alignments at radio wavelengths or the presence of resid¬ 
ual biases in the JVAS/CLASS 8.4-GHz radio polarization samples. 
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APPENDIX A: STATISTICAL TESTS 

To investigate the polarization RA correlations and to address the question of 
uniformity of their distributions, we use several stati s tical t ests in this paper. 
T hese are the Fourier metho d of lHawlev & Reebl^ il975|), the den s ity tes t 
of iRelgrims & Cud^ i2014b and the S and Z tests of lHutsemekersI il998h . 
For the sake of completeness and as they are little-known tests, we hereafter 
present a brief description of those. 


A1 The Hawley-Peebles Fourier method 

The Hawley-Reebles test taawlev & Reeble^ll975h is one of the com¬ 
mon methods used to study the alignments of galaxies. It is based on fit¬ 
ting the observed distribution of RA’s by a model of the form N(0i) = 
TV (1 + Ai cos 29i + A 2 sin 29i) where TV is the mean of the number of 
objects per bin and N(9i) is the observed number of objects in the bin 
centred in 9i. The number of bins is a free parameter. Ai and A 2 are the 
coefficients of the wave model which describe the degree of deviation of the 
distribution from being uniform. If the RA’s are not uniformly distributed the 
mean position angle is given by 0 = (1/2) arctan (A 2 /A 1 ). Agood mea¬ 
sure of departure from uniformity is the total amplitude A^ = Ai^H-A 2 ^. 
As easily understood, the bigger the value of A, the less uniform the distri¬ 
bution. The probability that the total amplitude exceeds by chance a given 
value of A is computed to be approximately Php = exp nA^/4^ 
where n is the number of objects in the sample. However, as far as small 
samples are considered, random simulations are required as the distribu¬ 
tion of A differs from a norm al Gaussian and a s this approximate relation 
fails far out in the wings fsee lGodtowskill? 012 l for a detailed discussion). 
For the simulated samples, polarization RA’s are uniformly generated and 
distributed among the sources. The probability Pup is then simply given 
by the percentage of random realizations having a A value higher or equal 
than that of the data. In this work, we tested all reported probabilities using 
random simulations and we only found marginal differences compared to 
those given by the approximate relation. These are actually smaller than the 
valuations caused by the choice of the number of bin. We decided to report 
only the probabilities computed through the Gaussian approximation. It is 
worth to mention that this statistical test is also dependent of the coordinate 
system in which RA’s are defined. However, as we use this test to study the 
uniformity of the polarization RA distributions inside relatively small re¬ 
gions of the sky and that the declinations of these regions are not too high, 
the changes are expected to be small. 


A2 The density test of Pelgrims & Cudell (2014) 

The intrinsically coord inate-invariant statistical test developed in 
iRelgrims & Cudefil j2Q14h introduces a polarization space (a 2-sphere) and 
treats polarization as points on it. The study of the density of these points 
allows a direct highlight of the direction for which an unexpected over¬ 
density is observed. For a given sample of sources a preferred direction 
is thus automatically allocated, to which corresponds a semi-analytically 
computed probability. This local probability, denoted by Pmin and being 
the p-value of a Roisson-Binomial distribution, gives the likelihood that the 
observed density toward the given direction is due to chance, i.e. assuming 
a uniform distribution for the polarization RA’s. A global probability, 
computed through Monte Carlo treatment, is also introduced and is denoted 
by . For a given sample, this global probability answers the question of 
the likelihood of a de tected significant over-den sity without concern of its 
direction. We refer to lRelgrims & Cud^ i2014h and lRelgrimj i2014l) for a 
complete description of this new test and for a discussion. Note however 
that throughout the study of the density of polarization points on the sphere, 
spherical caps having half-aperture angle 77 are used. This angle is the only 


parameter of this test and is strongly related to both the sky distribution of 
the sources and the strength of the alignment of the polarization vectors. 
It is therefore necessary to explore a large range of value of 77 . However, 
as this is the case with the bin width value for simple binomial tests on 
histograms, this method undergoes edge effects. A search for the veiy 
optimal value of 77 is therefore ambiguous. 

During our analysis (e.g. in Section [TH we report the mean position 
angle 9pQ for a given distribution. It is computed with respect to the spher¬ 
ical basis vectors *^he location given by the normalized 

vectorial sum of the source positions. This quantity has a meaning only 
when the maximum angular separation between studied sources is not lai'ge 
and when the position angle distribution is not uniform. 


A3 S and Z statistical tests 

A3.1 The S test 

The S test was developed by Hutsemekers (1998) in order to detect and 
statistically characterize the alignment features of polarization vector ori¬ 
entations which were first visually detected. This test is based on dispersion 
measures of FA’s for groups of riy neighbouring sources among the sample. 
For each object, the quantity 

riy 

d( 6 ») = 90-—V |90-|6)fc-6>||, (AI) 

riv 

k=i 

is computed, where the 9f.’s aiQ the polarization RA’s of the objects of the 
group, including the central one, in degree. For the object i, the mean dis¬ 
persion of the RA’s of the riy objects is computed to be the minimum value 
of d{9) and is denoted Si.lf n is the size of the whole sample under con¬ 
sideration, the statistic with the free parameter riy is defined as 

n 

Sd = — Si. (A2) 

n Z—^ 

i=l 

Sd measures the concentration of angles for groups of Uv objects close 
to each other in space (in 2 or 3 dimensions, as we shall see later). If the 
polarization vectors are locally aligned, the value of Sd will be smaller than 
in the case where the RA’s are distributed following a uniform distribution 
on the objects. 

The significance level (S.L.) to which one may assign the observed 
alignment to chance is evaluated through simulated datasets. Those are cre¬ 
ated by generating RA’s according to a uniform distribution or by shuffling 
the observed RA’s among the sources of the sample. This second procedure 
is found to be more appropriate to the detecti on of correlations b etween po¬ 
larization RA’s and the source locations fsee. lHutsemekerslll998h . This pro¬ 
cedure is chosen in this work each time we resort to dataset generation for 
the S and Z tests. The percentage of simulations with a value oi Sd lower 
than that of the data gives the S.L., i.e., the probability that the distribution 
of position angles is due to chance. 


A3.2 The Z test 

The Z test i s a non-parametr ic test originally introduced by Andrews & 
Wasserman iBietenh(^l 1 98d) to quantify the correlation between the RA’s 
and the position of sources on the sky. Its basic idea is to compute for each 
object i, the mean direction 9i of its riy neighbours, excluding this time the 
central object i, and to compare this local average to the actual polaiization 
RA of the object i, namely, 9i. Specifically, the RA’s of the Uy nearest neigh¬ 
bours around each object i but excluding the latter, are used to computed the 
mean resultant vector 

, (A3) 

where the factor 2 accounts for the axial nature of the polarization. The 
mean direction 9i is given by the normalized mean vector Y i through 


/ Tly riy 

Yi = — ( ^cos20fe, ^ 
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yi = ^cos20i, sin20i) . (A4) 

A measure of the closeness of 6i to 9j is the inner product Dij = • Yj, 
where = (cos29i, sin20i). If the PA’s are con'elated to the source 
positions, then, on average, 9i will be closer to 9j=i than to 9j^i which, in 
turn, implies Dij to be larger for i = j than for i ^ j. An approximately 
nor mally distributed statistic Zc is computed from these quantities as follow 
(cf. iBietenhofll 19861) : 

n 

Zc=-\^ Zi (A5) 

n 

i=l 


where 

^ _ ri - (n + 1) /2 

— -/ - 

VW12 


(A6) 


and where ri is the rank of Dij=i, when the Dij=i^ri’s are sorted in 
increasing order and n is the size of the studied sample. As correlations of 
PA’s with source positions induce large Dij=i, Zc is expected to be larger 
than zero if polarization vectors are coherently aligned. 

Despite the fact that the Zc statistic is approximately normally dis¬ 
tributed, simulations are needed to evaluate the S.L. because of the mutual 
dependence of the Dij=i^n, especially for large values of ny. Therefore, 
as for the S test, the S.L. is given by the percentage of simulations that show 
greater value of Zc than the one corresponding to the observations. 

_ A modification of the Andrews & Wasserman test was proposed by 

iHutsemekerj il998h . Instead of considering the inner product with Yj, he 
proposed to compute Dij = • Yj. This definition gives more weight to 

the groups of sources having similar PA values. Indeed, aligned polarization 
vectors imply a large norm of j which provides a natural measure of the 
dispersion of the position angles and leads to a large Dij. Apart from this 
variation, the modified statistic is computed in the same way as the original 
test of Andrews & Wasserman. In this work, we use this modified version 
when we refer to the Z test as it is more sensitive to local alignments and 
thus, more adapted to the search for such features. 
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ASA One free parameter statistics and physical interpretation 

In both S and Z tests, the number of nearest neighbours riv is a free parame¬ 
ter which has to be explored. Indeed, this parameter is not devoid of physical 
meaning as it is related to a characteristic scale of the nearest neighbouring 
groups, in two or three dimensions. This would be strictly true for a sample 
uniformly scattered across the whole sky. As the observed samples show 
deviations from homogeneity and, more importantly, as only a pai't of the 
entire celestial sphere is covered, this parameter does not show a straight 
correspondence with a typical size of groups. A dispersion is naturally ex¬ 
pected. Nevertheless, if correlations between polarization orientations occur 
for a typical scale or if some sub-samples well delimited -in space- present 
such alignments, it is clear that the S.L. will be smaller for the correspond¬ 
ing value of Uv than for others. Therefore, to test the uniformity of the 
polarization orientations and explore their characteristics, it is necessary to 
estimate the S.L. across a wide range of values of riy. To build the groups 
of nearest neighbouring objects, the relative distances between sources is 
required. For this purpose, and assuming a flat and isotropic Univer se, we 
consi der the line-of-sight comoving distan ce r to the observer (e.g. iHogd 
I 1999 I) . In agreement with lAde et alJ j20l4) . our calculation are made with 
Hq = 68 kms“^Mpc“^ and = 0.31. 

The relative distances between sources are then simply computed 
through the use of their rectangular coordinates and the groups of near¬ 
est neighbours constructed. A three dimensional analysis is obviously only 
feasible for samples for which redshift measurements are available. Nev¬ 
ertheless, a two dimensional one, considering nearest neighbours on the 
celestial sphere rather than in the 3D space, is not devoid of interest and is 
applicable for all samples imposing r = 1. 


AS.3 Parallel transport for coordinate-invariant statistics 

The PA’s of quantities which are transverse to the line-of-sight to their corre¬ 
sponding sources, i.e. projected in the planes orthogonal to these directions, 
are dependent of the coordinate system in which the source positions are re¬ 
ported. As the changes of PA values depend on the source positions, results 
of statistical tests change from coordinate syste m to coordinate system. In 
order to overcome this coordinate dependence, Ijain et all l2004h modified 
the S and Z tests introduced above. 

Instead of computing statistics directly from the PA’s each evaluated 
with respect to their own meridian, they introduced correction to the PA’s 
that involves the relative sky positions of the sources. In order to compare 
in a coordinate- invariant manner the PA’s of two sources separated on the 
celestial sphere. Uain et aljj2004h did parallel transport the polarization vec¬ 
tor of one source to the other. The parallel transport is performed along the 
sphere geodesic passing through the positions of the two sources. The idea 
is thus to add a correction term to the PA values, this correction being sim¬ 
ply the difference between the angles that forms the geodesic with one of 
the basis vectors at the locations of the two sources. Therefore, while eval¬ 
uating Eg. lAll or Ea. lA3l for each central object i, the polarization PA’s of 
the Tiy neighbouring sources are transformed as 

9'f^=9k + , (A7) 

where each A^-^i is the correction due to the parallel transport of basis 
vectors from th e position of the source k to the position of the source i (see 
Ijain et alJ20Q4l for further details). 

Apart from this modification, the statistical tests are performed as ex¬ 
plained abov43 arid lead to coordinate-invariant results, as wanted. 

® Note that I Jain et alJ j2Q04h used another definition for the statistics Sd 
when evaluating Eq. IaTI 
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